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We demonstrate a convenient approach for simultaneously manipulating the amplitude and polarization of light beams by 
means of the modulation of the correlation structure. As an illustration, we constructed a periodic correlation structure that 
can generate an arbitrary radially polarized array (RPA) beam of a radial or rectangular symmetry array in the focal plane 
from a radially polarized (RP) beam. The physical realizability conditions for such source and the far-field beam condition 
are derived. It is illustrated that the beamlet shape and the state of polarization (SOP) can be effectively controlled by the 
initial correlation structure and the coherence width. Furthermore, by designing the source correlation structure, a tunable 
OK-shaped RPA beam and an optical cage are demonstrated, which can find widespread applications in non-destructive 
manipulation of particles and living biological cells. The arbitrariness in the design of correlation structure prompted us to 
find more convenient approaches for controlling the statistics of light beams in terms of amplitude and polarization. @ 2016 


AIP Publishing LLC. 


In the past decade, array beams have seen a rapid growth of interest in the use of optical interconnections, ultracold 
atoms and optical tweezers [1—4]. Therefore, a variety of approaches have been introduced for generating array beams such 
as arrays of microlenses, binary phase gratings, and interferometric technique [5-7]. In general, array beams in previous 
studies including Gaussian array beams, dark-hollow array beams, flat-topped array beams, and Bessel beams with a spatially 
homogeneous state of polarization (SOP), are within the framework of scalar beams [8—10]. However, if the SOP of a light 
beam is inhomogeneous but is instead a spatially varying polarization distribution in the transverse beam cross-section, 
vector theory is needed. Due to their essential features that are not present in scalar beams, vector beams have started to open 


up new research opportunities in optics [11]. 
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On the other hand, significant efforts have been devoted to the generation of cylindrical vector (CV) beams, a family of 
vector solutions of Maxwell’s equation, which features a cylindrical SOP orthogonal to the propagation direction. As an 
important class of CV beams, radially polarized (RP) beams have been extensively investigated due to their special features 
such as a tight focal beam spot and strong longitudinal component of the electric field at the focal plane, which have proved 
to be a powerful tool for many applications in polarization information encryption, optical data storage, super-resolution 
imaging, material processing, plasmonic focusing and optical tweezers [12—23]. Various methods for generating RP beams in 
terms of interferometric technique, modified laser cavities, binary optical elements, half-wave plates, and liquid crystal phase 
modulators have been widely explored [24-32]. However, very few work has referred to the generation of radially polarized 
array (RPA) beams. In some cases of practical application, for instance, multiple manipulation of nanoparticles, living 
biological cells and organelles within cells, a RPA beam have significant advantages in comparison with conventional scalar 


beams [3, 12]. 


The notion of coherence is an important fundamental of modern optics. It is well known that the coherence is 
intrinsically connected with the spectral intensity, SOP, beam quality and the orbital angular moment (OAM) [33-42]. For 
instance, one cannot determine coherence from the visibility of interference pattern when the SOP of a vector beam is 
ignored. Recently, a sufficient condition was introduced by Gori et al. for devising genuine spatial correlation functions on 
taking into account the non-negative definiteness constraint of the cross-spectral density (CSD) [43, 44]. A variety of light 
beams with peculiar correlation properties have been proposed, for which unexpected features have been revealed such as 
self-accelerating, ring-shaped profile as well as optical lattices in the far field [45-56]. In this paper, we extend the work from 
scalar case to the vector case, and demonstrate the ability of correlation structure to simultaneously manipulate the amplitude 
and polarization. By constructing a tunable periodic correlation structure, we have achieved a physically feasible approach 
for converting a RP beam into an arbitrary RPA of highly directional beam in the focal plane. The realizability conditions for 
such source and the beam condition for radiation generated by such source are established. It has been shown that the far- 
field SOP and the beam shape can be effectively manipulated by taking advantage of the source correlation structure and the 
coherence width. As further explorations, a tunable OK-shaped RPA beam and an optical cage are demonstrated by using a 
similar approach, which will be found important applications in non-destructive manipulation of living biological cells and 
capture of metal nanoparticles. Owning to the rapidly growing current interest in an unified theory of coherence and 
polarization, it is hoped that the arbitrariness in the design of correlation structure will stimulate further research into the 


generation and manipulation of novel vector beams. 


In general, the spatial coherence properties of a light field can be characterized by means of CSD. Within the paraxial 
approximation yield, the second-order correlation properties of a quasi-monochromatic, statistically stationary 
electromagnetic beam propagating along z-direction can be described by the 2x2 CSD matrix in space-frequency domain, 


whose elements are given by [33] 


Vag (rr) = (E; (1,)E;(r,)); (a,B=xy), D 


where r= ( X, y) (j=1,2) is the arbitrary position vector in the source under the Cartesian coordinate. The angle brackets 
denote time average and the asterisk represents the complex conjugate. E (æ =x,y) denotes the fluctuating electric field 


component along a-direction. For simplicity, we not specify the frequency @ throughout this paper. 


It is well known that the coherence is closely associated with the polarization. Naturally, one can expect to modulate the 
polarization by means of the coherence. However, although the correlation model of light beam has rich variety of choices, it 
cannot be chosen at will because of the non-negative definiteness constraints [43]. To ensure that the constructed correlation 
structure is reasonable and physically realizable, the elements of CSD matrix have to be non-negative definite kernels. Let’s 


define the quadratic form 


O=S [AMAT Gm dnd’, (a, B=xy), a 
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where f, is an arbitrary (well-behaving) function. Then, the non-negative constraint means that, for any choice of f and f, 
Q must be non-negative. However, to make sure that Q meets the conditions for a given form of Ta is quite difficult, 
because of the arbitrariness in the choice of f and fy A sufficiency condition for satisfying the non-negative constraint in Eq. 


(2) has recently been proposed [43]. It is shown that a well-defined form for the elements of CSD with the help of the 


reproducing kernel Hilbert spaces can be expressed as an integral 


Dog (t0y) = Í | Pp) KUG,. VK pm, vdv’, (a, 8 =x, y), (3) 
where V=(V,5V,)> Pag is an non-negative weighting functions and can consist of a set of delta functions with positive 
coefficients. Ķ is an arbitrary kernel and can be regard as the physically available response functions associated with the 
electric field components x. On inserting from Eq. (3) into Eq. (2), for any non-negative definiteness constraint of Q, 


Pag should satisfy the following inequality [44] 


Px (y) > 0, py (Y) 0, and p,.(v)P,,(V) -|p 0, (4) 


for any v. As a consequence, the elements of CSD matrix constructed by means of the recipe given by Eq. (3) are necessarily 


non-negative definite, namely, it is physically realizable, as far as Eq. (4) is satisfied. 


Since K represent arbitrary kernels of linear transformation associated with v, it has various choices such as Fourier, 


Laplace, Hankel, mellin, and etc.,. Here, for simplicity, let us suppose kernel K, to be a conventional Fourier-like form, 


K (Œ, v) =U,,(r)exp[-2zig,(r)- v], (5) 


where Ọ „(r)is possible complex profile function and g, is arbitrary vectorial real function. On inserting Eq. (5) into Eq. (3), 


we get the following form 


TTT) =U GU E, ) Pap | 850) — 8.0) |- (6) 


where the tilde represents the Fourier transform. It is seen that p 5 coincide with the source complex degree of coherence 


(DOC). Since g, can be chosen as will, the elements [.,,can have rather sophisticated structures. Let us consider the 


ap 
simplest case, for instance g (r) =r. In comparison with Eq. (3) and Eq. (6), one finds that the source correlation function is 


the Fourier transform of weighting function p, g 
It is worthwhile to recall a typical case that the correlation structure has a Gaussian form (known as Schell-model (SM)) 
[33—45]. Then, the weighting function p, 7 has the following form 


(SM) 


Pop (V) = 2705 exp| —27’o?v? |. (7) 


In the following, let us construct p, pto have a M x N radial array structure with the help of a periodic comb function, 
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where C,=M (2N +1) and 2, = 2an/ (N + 1) sT is the radial of the inner ring. On the other hand, since profile function can be 


chosen at will, let us take into account U(r) having a RP form [12, 13], 


ULU r) = Že ol 1 t) (a,B=x,y), (9) 
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where w denotes the transverse beam size. On inserting from Eqs. (8) and (9) into Eq. (6), we obtain the following 


expressions for the elements of CSD matrix 
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One can notice that Eq. (10) characterizes a RP beam with a Schell-model array correlated (SMAC) function of radial 


symmetry. By evaluating Eq. (8), it is easy to find 


Pa (Y) = Pp (Y) = p, (Y) > 0, and |p,,(v)} = p,.(v)p,,(¥). (11) 


According to the Eq. (10), it is clear that Eq. (4) is satisfied for any v. Furthermore, for this new type of non-uniformly 


correlated RP beam, the beam condition for radiation generated by such source can refer to Refs. [49—52], and it is derived as 
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FIG. 1. Contour graphs of the source spectral DOCs and the corresponding cross lines for radial and rectangular SMAC RP beams as a 
function of g, . (al)-(b1): o, =Imm,2a =r, =1mm, (@2)-(b2): o, =10mm,a = 2r, = 1mm. 


Figure 1 shows the source spectral DOCs and the corresponding cross lines for radial and rectangular SMAC RP beams 


as a function of >: The parameters are chosen to be A= 632.8nm, w, = Imm, M =2,N=7,M'=N'=1, and o, represents the 


initial coherence width. It can be seen that the source correlation array structures are determined by the periodic function in 
the spectral DOC, and the region of correlation array is controllable owing to the Gaussian correlation function. As the 


increase of the value g, the correlation array region gradually expand. Owning to the nonconventional correlation structure, 


light beams thus exhibit very different modulation characteristics. We shall emphasize that the propagation-induced statistics 


changes are completely due to the correlation structure and do not result from any deterministic phase modulation. 


In the following, let us consider the paraxial propagation of this beam through a stigmatic ABCD. The components of 


spectral intensity at any point p = (x, y) across the transverse plane can be evaluated by the generalized Collins formula [13, 33] 


ikA ik 
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where / is the wavelength and k =277/A denotes the wave number. On substituting from Eq. (10) into Eq. (13), we obtain 


the following expressions for the components of spectral intensity 
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By applying a coordinate transform, we further obtain the following components of spectral intensity for a (2M' +1) 


x(2N' +1) rectangular SMAC RP beam, 
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with 


6, =x—AmB/2a, €,=y—AnB/2a, and C, =(2M'+1)(2N' +1). (19) 


where ais the separation distance. It is easy to find that the realizability conditions for a rectangular SMAC RP source is 


coincide with the radial SMAC RP source, and the beam condition is obtained as follows 


K 2 
E a Í a -max(M.N)] (20) 
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Let us now take into account the focusing vector properties. Here, we assume that the source beam propagates a distance 


of/=150mm, and then focused by a thin lens with focal length f =150mm. The elements of the transfer matrix between the 


source plane and the focal plane are given by 


A=1-~, B=f, Cas, D=0. 


(21) 
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FIG. 2. Contour plots of the SOPs in the source plane for radial and rectangular SMAC RP beams as well as that in the focal plane as 
function of g, . (a2)-(b2): o =1mm, 2a = =1mm, (a3)-(63): o =10mm,a = 27, =1mm. 


Figure 2 shows the contour plots of the spectral SOP for radial and rectangular SMAC RP beams in the source plane as 
well as that in the focal plane. It is clearly seen that the spectral density in the focal plane exhibits array distribution with 


radial or rectangular symmetry. This means that the Fourier-like relation between the initial correlation function and the far- 
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field spectral intensity distribution is still valid for vector beams [56]. At the same time, each beamlet evolves from a 


Gaussian shape into a dark-hollow profile when the value of, increases. In addition, each beamlet copies the initial SOP and 


displays RP structure. This implies the initial overall polarization structure is completely broken induced by the 
nonconventional correlation structure. More explicitly, polarization modulation and amplitude modulation are not 
independent features that can be simultaneously manipulated by means of the specially designed correlation structure. Our 
further study shows that this reproduction induced by the source correlation structure is also applicable to other vector beams 


such as generic CV beams. However, the present example for modulating a vector beam is more than just a simple split (as it 


might appear) of initial SOP, and needs for further investigation. 
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FIG. 3. Contour plots of the SOP, corresponding intensity components and the cross lines of a radial SMAC RP beam in the focal plane 
with different Oy: (al)-(a3): o =1mm, (b1)-(b3): o, =10mm. 


Figure 3 illustrates the SOP and the corresponding intensity components of a radial SMAC RP beam in the focal plane as 


a function ofg,. Previous studies show that a separable weighting function p, P in the Cartesian coordinate system leads to a 


rectangle-like far-field intensity profile [45, 56]. However, Fig. 3 shows that this is not a sufficient condition and a circular 
beam can be obtained as well, even p, pis derived under the Cartesian coordinate system. Different from the polarization 
reproduction shown in Fig. 2, it can be seen that the SOP in Fig. 3 displays a more complicated structure. In the case of lower 


value of coherence width such aso, =1mm, one finds from Figs. 3(a1)—3(a3) that the most region of the beam is azimuthally 


polarized, except the edge portion displays a RP structure. However, as the increase of the coherence width such 


aso, =10mm, Fig. 3(b1) shows that the SOP in the central region evolves into RP structure. Furthermore, the corresponding 


intensity components, as presented in Figs. 3(a2), 3(a3), 3(b2) and 3(b3) clearly demonstrate this variation in polarization 


state. 


To demonstrate more about this method, Fig. 4 illustrates an OK-shaped RPA beam in the focal plane based on the 
special designed initial correlation structures. It is interesting to note that not only the array structure can be arbitrarily 
designed, but also the SOP can be precisely controlled by means of the initial correlation structure. The ability of generation 
arbitrary vector array beam can find widespread applications in particle manipulation and optical engineering. By using a 
similar approach, Fig. 5 demonstrate an optical cage intensity profile around the focal region. In particular, the size of optical 
cage can be effectively regulated by taking advantage of a thin lens. Such an optical cage is essential needed for non- 
destructive manipulation of biological cells and capture of metal nanoparticles [2, 3]. Finally, it is important to point out that 
the correlation structures mentioned above can be generated experimentally by means of a liquid crystal spatial light 


modulator or a method proposed recently in Refs. 44 and 45. 
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FIG. 4. Contour plots of an OK-shaped RPA beam and its corresponding intensity components in the focal plane. 
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FIG. 5. An optical cage is derived around the focal region. 


In conclusion, we have demonstrated a convenient approach for simultaneously manipulating the amplitude and 
polarization of light beams by means of the design of the correlation structure. As an illustration, by constructing a periodic 
correlation structure, a RP beam is focused to an arbitrary RPA in the focal plane. The realizability conditions for the 
correlation structures and the far-field beam condition are established. It has been shown that the far-field spectral intensity 
distribution and the SOP can be effectively controlled by modulating the initial coherence width. Furthermore, correlation 
structures give rise to a tunable OK-shaped RPA beam and an optical cage are demonstrated. The ability to appropriately 
manipulate amplitude and polarization of light beams by designing the initial correlation structure gives an incentive for the 
generation of more complex vector beams such as high-order vector beams, generic CV beams or Airy beams and promises 


important supports in the fields of imaging, laser machining, and optical tweezers. 
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